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RESULTS FROM FLIGHT AND SIMULATOR STUDIES OF A 


MACH 3 CRUISE LONGITUDINAL AUTOPILOT 


Glenn B . Gilyard and John W . Smith 
Dryden Flight Research Center 


INTRODUCTION 


For subsonic jet transports, accurate control of altitude and Mach number is 
necessary for such items as air traffic control and good ride qualities . For aircraft 
operating at high altitudes and supersonic Mach numbers , accurate control becomes 
increasingly necessary for maximum range performance. The aircraft of the YF-12 
series present an excellent opportunity to document the interaction of various auto- 
pilot schemes , upper atmospheric effects , and aircraft dynamics for the advancement 
of supersonic cruise aircraft technology . 

The precise control of flightpath and speed becomes increasingly difficult at 
high altitude , high speed flight conditions . Large excursions in altitude or speed , 
or both, have been experienced with the XB-70 (ref. 1) , YF-12 (ref. 2) , and 
Concorde (ref. 3) airplanes. Factors that contribute to such performance include 
decreases in aircraft stability, low static pressure, atmospheric disturbances such 
as temperature variations , and unfavorable aircraft-propulsion system interactions . 
Low atmospheric static pressure makes precise altitude control difficult because of 
poor pressure transducer resolution . In addition , the combination of high altitude 
and high speed contributes to an unfavorable balance between kinetic and potential 
energy , which in turn makes large altitude changes necessary to correct for small 
Mach number changes when the aircraft is flown in a Mach hold mode through the 
elevator . Interactions between the aircraft and the propulsion system also become a 
significant factor at high altitude , high speed flight conditions . Reference 2 presents 
flight data illustrating that longitudinally, the YF-12 aircraft's phugoid mode has a 
larger amplitude and less stability with the inlets operating automatically than with 
the inlet geometry fixed . 

One aspect of the YF-12 research program is to identify the deficiencies of the 
original autopilot, which is approximately 10 years old, and to implement an auto- 
pilot improvement program. The first part of the improvement program was to update 
the existing altitude hold mode of the autopilot with minimal modification to reflect 
recent knowledge and experience. In subsequent phases, an autothrottle for speed 
control is to be implemented and an integrated airplane-propulsion system controller 



(ref. 2) is to be developed to utilize previously identified lateral-directional 
airplane-propulsion system interactions (refs. 4 and 5). 

This report presents flight data obtained with the original altitude hold and Mach 
hold control systems and the results of simulator studies , which include parametric 
variations of various gains . Flight test results for the modified altitude hold autopilot 
and the Mach hold autopilot when connected to an angle-of-attack-insensitive static 
pressure source are also presented . A condensed version of this report is presented 
in reference 6 . 


SYMBOLS AND ABBREVIATIONS 


Physical quantities in this report are given in the International System of Units 
(SI) and parenthetically in U.S. Customary Units. The measurements were taken in 

Customary Units. The decibels in this report are 20 ’ 


A amplitude of sine wave oscillation 

A(N) actuator nonlinearity , ±0. 15° hysteresis 

Aj(s) inboard actuator transfer function 

Aq (s) outboard actuator transfer function 


ADC 


B 



D 

DI 


F 

Fhp(s) 


air data computer 

normal acceleration at center of gravity (unless otherwise noted) , g 

first bending mode 

first bending mode transfer function 

nonlinearity width 

X2-X1 

damping index, — y ■ (derived in app. B) 

Xq X^ 

AM/M 

high-pass filter 
structural filter 
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F (s) .F 2 (s) shaping filter 


D 


GyCN) 


G'(s) 

H(s) 

H'is) 

h 

i 

K 

M 

m 

PLA 


PgCa) 


si 


frequency, Hz 
describing function 

backup gyro nonlinearity 

equivalent system forward-loop elements 
SAS feedback loop filters 
equivalent system feedback elements 
altitude, m (ft) 

gain 

Mach number 

peak number of sine wave oscillation 
power lever angle, deg 

2 2 

static pressure, N/m (Ib/ft ) 


variation of static pressure error with angle of attack, N/m'^/deg 
(Ib/ft^/deg) 

2 2 

sea level static pressure, N/m (Ib/ft ) 


total pressure, N/m^ (Ib/ft^) 


R 

SAS 

s 


2 2 

impact pressure , p. ~ P„. N/m (Ib/ft ) 

'2 ® 

PtJPs 


stability augmentation system 
Laplace operator , 1/sec 
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30 


X(t) 

X 

m 


6 

6 

e 

? 

0 
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^(s) 
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X 

<P 

CO 

CO^ 


(O 


n 


time, sec 

time required to reach and maintain altitude within ±9.1 meters 
(±30 feet) of steady state, sec 

= Ae " 

peak amplitude of sine wave oscillation; m = 0 is the reference peak 
and m = 1 , 2 , 3 , . . . are the successive peaks 

wing reference plane angle of attack , deg 

= -c/Vi - 

average eleven deflection , deg and rad 
autopilot gains 

damping ratio 

pitch attitude , deg and rad 

controlled attitude, rad 

reference attitude, rad 
aircraft transfer function 
time constant 

time constant referenced to sea level 
autopilot time constant 

phase angle , deg 
frequency, rad/ sec 
damped frequency, rad/sec 

natural frequency, rad/sec 
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Subscripts: 

1,2 compensated static pressure location . 

3 uncompensated static pressure location 

A dot over a quantity denotes the time derivative of that quantity , 


AIRCRAFT DESCRIPTION 


The YF-12 airplane is an advanced, twin-engined, delta- winged interceptor 
designed for long-range cruise at Mach numbers greater than 3.0 and altitudes 
above 24,400 meters (80,000 feet) . A photograph and a three-view drawing of the 
airplane are shown in figures 1 and 2 . Pertinent physical characteristics are given 
in reference 4. 



E-23090 

Figure 1. YF-12 airplane. 



Figure 2. Three-view drawing of YF-12 airplane 



Two nacelle-mounted, all-movable vertical tails provide directional stability and 
control . Each vertical tail is canted inward and pivots on a small stub section 
attached directly to the top of the nacelle . Three ventral fins enhance directional 
stability; one is fastened along the fuselage centerline and two small fixed fins 
extend down from the nacelles below the two main vertical tails . The centerline 
ventral fin folds out of the way for landing . 

Two elevons on each wing , one inboard and one outboard of each nacelle , perform 
the combined functions of ailerons and elevators . 

The airplane has two axisymmetric , variable-geometry, mixed-compression 
inlets, which supply air to two J58 engines. Each inlet has a translating spike and 
forward bypass doors to control the position of the normal shock in the inlet . An 
automatic inlet control system varies the spike and bypass door position to keep the 
normal shock in the optimum position . Manual control of the spike and bypass doors 
is also available to the pilot . 


ELEMENTS OF FLIGHT CONTROL SYSTEM 


The main elements of the pitch-axis flight control system are the air data sensors , 
the air data computer , and the original autopilot control system . These items are 
briefly discussed in this section; a more detailed description is presented in 
appendix A . 


Air Data Sensors 

The air data parameters , static and total pressure , are obtained from the compen- 
sated nose boom illustrated in figure 3 . A more detailed view of the three static 
pressure source locations is presented in figure 4 . The air data computer is con- 
nected to p , which is compensated to minimize position error corrections , but 

is sensitive to angle of attack as shown in figure 5 . 

The variation of static pressure error with angle of attack in figure 5 , which is 

referred to throughout this report as the nominal p (a) , was determined from slow 

s 

pullup-pushover maneuvers of approximately 20 to 30 seconds' duration. The 

nominal p (a) at Mach 3.0 and 23,600 meters (77,500 feet) at a typical trim angle of 
s 

attack is equivalent to an altitude error of 49 meters (161 feet) per degree angle of 
attack or a Mach error of 0.011 per degree angle of attack. 

The p source was not originally on the nose boom , but was added for research 
®3 

purposes . The p measurements require large position error corrections tran- 
®3 

sonically but have negligible sensitivity to angle of attack . 
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Hemispherical head sensor details 


Figure 3 . Three-view drawing of nose boom and pitot-static probe showing 
hemispherical head flow direction sensor. Dimensions in centimeters (inches). 
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Figure 4. Static pressure source locations on compensated nose boom. 
Dimensions in centimeters (inches). 
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Figure 5. Variation of static pressure error with angle of attack, 

p (a). M^3.0, h = 23, 600 meters (77 , 500 feet) . At the trim point, 
s 

2 2 

Ap /Aa = -23.9 N/m /deg (-0.50 lb/ ft /deg) , which is equivalent 
to Ah/Aa = 49 m/deg (161 ft/deg) or AM/Aa - 0. 011/deg. 


At 23,600 meters (77,500 feet) , the approximate lags due to the tubing lengths 
of the static pressure and total pressure systems are 2.00 seconds and 0.75 second, 
respectively . 


Air Data Computer 

The air data computer (ADC) is an electromechanical device that receives total 
and static pressure from the nose boom and computes Mach number and altitude 
information for the cockpit display and the autopilot. The threshold of the static 
pressure loop at an altitude of 23,600 meters (77,500 feet) is approximately 3.7 meters 
(12 feet) and its frequency response is constant to approximately 0.5 radian per 
second . 


Original Altitude and Mach Hold Control Systems 

Figure 6 is a pitch-axis block diagram of the vehicle control system. At Mach 3.0 
and an altitude of 23,600 meters (77,500 feet), the gains and time constant for the 
original autopilot, which are referred to throughout this report as the nominal 
settings, are as follows: 


Kq = 440 deg Q/ln p 
efi s 
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I 



I 


Xq = 170 deg 9/Aln 
Xq = 14 (deg/sec 0)/AZn p„ 

Vh ® 

X„ = 20 (deg/sec §)/Am (R - 1) 
^fM 

X = 61 deg e/AIn (R - 1) 


Xp =0 deg 9/ln (R - 1) 
^Af 


^0 


= 1 


6a = 0.295 deg 6 /(deg/sec 6) 

tJ 6 


t 


5.75 


6p = 3 . 50 deg 8 /deg 0 
fc»L e 

The pitch stability augmentation system (SAS) is used full time for normal aircraft 
operation and attitude hold is the primary autopilot mode . In addition to the basic 
attitude hold autopilot , the pilot can select either altitude hold or Mach hold as an 
outer loop of attitude hold . 

The altitude hold autopilot receives altitude rate and altitude error information 
from the ADC and commands attitude changes proportional to the altitude rate, 
altitude error , and integral of altitude error . Similarly , the Mach hold autopilot 
receives Mach error information from the ADC and commands attitude changes 
proportional to the Mach error and the integral of Mach error . 

Both the pitch SAS and attitude hold loops, which are inner loops for the altitude 
hold and Mach hold modes , performed satisfactorily . 


FLIGHT INSTRUMENTATION 


A standard set of stability and control parameters was recorded. The angular 
rate and linear acceleration instrumentation was alined with the body axes . The 
angle-of- attack indicator was attached as a dogleg on the nose boom , which also 
contained an airspeed-altitude probe (fig. 3) . A fixed four-port pressure- sensing 
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Schedules are given in figure 52. 


hemispherical head was used to obtain 
angle-of-attack measurements on the 
airplane (ref. 7) . The lag associated 
with the hemispherical head was signifi- 
cantly different from that predicted by 
first-order lag theory (ref. 8). The 
actual angle-of-attack system lag at 
Mach 3 . 0 flight conditions was approxi- 
mately 0.4 second. 

In addition , data were recorded at 
test points within the autopilot control 
system to analyze system performance 
and to detect system problems . All the 
data were recorded on magnetic tape at 
200 samples per second. 


SIMULATION SYSTEM 


A digital simulation was implemented 
to investigate closed-loop aircraft-control 
system problems . The simulation included 
aircraft and propulsion dynamics , the 
control system , and an air data computer 
model. An integration interval of 20 milli- 
seconds was used for the entire simulation. 


Aerodynamic and Propul sion Dynamics 

The simulation was modeled for the 
flight conditions at Mach 3 . 0 and an alti- 
tude of 23,600 meters (77,500 feet) . The 
simulation was a modification of that de- 
scribed in reference 9 and included the 
three longitudinal degrees of freedom . 
Dynamic pressure flexibility corrections 
were included for pitch control , although 
no structural modes were simulated. The 
variation of density with altitude was also 
included . 

High speed inlet operation and the 
afterburner range of engine operation were 
represented in sufficient detail to permit 
investigation of the aerodynamic and 
propulsion system interactions. Only 
normal inlet operation was modeled . 
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Control System 


With the exception of the speed stability system , which has no effect at the flight 
condition of the study , the control system illustrated in figure 6 was implemented . 
The schedules were programed as functions of static or dynamic pressure , or both . 
All control system dynamics above 5 hertz were not modeled. 


Air Data Comput er 

Accurate modeling of the ADC was required because the ADC generates the alti- 
tude and Mach signals input to the autopilot . Initial attempts at linearization proved 
inadequate in the complete closed-loop simulation . Figure 7 presents the complete 


In Pj at autopilot 
engagement 



Figure 7. Simulation mechanization of air data computer. 
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ADC model (based on the ADC description in appendix A) used in the simulation . 
The system characteristics were defined as follows: 

T„ =1.8 sec 
Pc 


Pc 


?p =1.0 
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s 


Ps 


X 


p 


= 0.2 sec 


“n =0-83V% 


? = 0.63 


K 


q 


c 



2 


q 


c 


The limits in the servo loops were due to the rate limits of the motor used in the 
rebalancing hardware . The slip clutch limits were due to the limited range of the 
potentiometer , which was designed to get maximum resolution at the expense of a 
larger error output range . 

The inputs to the ADC model were true static pressure , p , and total pressure , 

s 

p. . Angle-of-attack sensitivity was added to p and first-order line dynamics were 

h ® 

added to p and p. . 

« *2 


The outputs of the ADC were logarithmic functions , which were approximately 
linearly proportional to altitude and Mach number over the perturbation range of 
interest . 
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RESULTS AND DISCUSSION 


Altitude Hold Aut opilot 

This section presents the results of the altitude hold autopilot investigation . 
The discussion covers the original altitude hold performance; simulation studies, 
which include parametric variations of the feedback gains; the actual hardware 
changes made to improve the autopilot; and the in-flight performance evaluation 
of the modified altitude hold mode . 


Original Altitude Hold Operation 

The altitude hold mode of the YF-12 autopilot was designed as a loiter mode for 
use at altitudes of approximately 9100 meters (30,000 feet) , although operation was 
permitted at altitudes up to 18,300 meters (60,000 feet) . The control system does 
not limit altitude hold operation to 18,300 meters (60,000 feet) since the autopilot 
gains extend to the altitude and speed capabilities of the aircraft . However , the 
U.S. Air Force performed a preliminary evaluation of the altitude hold mode at 
altitudes exceeding 21,300 meters (70,000 feet) and found its operation to vary 
considerably from day to day . Occasionally , altitude could be held within ±15 meters 
(±50 feet) while at other times altitude diverged in an unacceptable manner. 

An example of acceptable altitude control at approximately Mach 3 . 0 and an 
altitude of 23,600 meters (77,500 feet) is presented in figure 8. The pilot described 
the atmosphere as stable as evidenced by the ease of maintaining Mach number and 
altitude conditions . The low frequency limit cycle , which has a period of approxi- 
mately 32 seconds, is due to the air data computer's static pressure threshold; 
the high frequency oscillatory characteristics, which can be seen in 6^, are due to 

the angle-of-attack sensitivity of the static pressure source, p . The normal accel- 

®2 

eration parameter at the center of gravity was noisy and sensitive to engine power 
level . 

An example of unacceptable altitude hold control at approximately Mach 3 . 0 
and an altitude of 23,600 meters (77,500 feet) is presented in figure 9. The pilot 
described the atmosphere as unstable because of the difficulty in maintaining alti- 
tude and the sudden increases in the indicated Mach number , On this occasion the 
atmospheric disturbances apparently induced the erratic altitude behavior; however, 
other types of disturbances or untrimmed conditions can initiate the autopilot insta- 
bility . The high frequency , divergent-convergent oscillations in this example are 
due to the angle-of-attack sensitivity of the static pressure ports, p (fig. 5) . As 

®2 

angle of attack increases , Ap /Aa becomes more negative and the high frequency 

®2 

motion diverges; as angle of attack decreases, Ap /Aa becomes less negative and 

®2 

the high frequency motion converges . 
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Autopilot 



0 12 3 4 

Time, min 


Figure 8. Good original altitude hold; stable atmosphere . 

In both the previous examples , the static pressure sensitivity to angle of attack 
is evident. In the latter example, the divergent, high frequency, large amplitude 
oscillations were disturbing to the pilot . 

A short evaluation of the altitude hold autopilot was conducted with the ADC 
connected to p , which is not sensitive to angle of attack . A typical example of 
®3 

altitude hold in this configuration is shown in figure 10 . The high frequency oscil- 
lations were eliminated as expected; however, the low frequency mode was diver- 
gent . The time history shows that after 2 1/2 minutes the autopilot reached satura- 
tion (±2.3® elevator authority) , as evidenced by the ramp type of elevon input 
caused by the autotrim feature of the autopilot . 

These examples demonstrate two major problems with the original altitude hold 
autopilot: the neutral to divergent damping of the low frequency mode and the high 
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Autopilot 



Time, min 


Figure 9. Unacceptable original altitude hold; 
unstable atmosphere . 

frequency , divergent-convergent oscillations due to the sensitivity of static 
pressure to angle of attack . 


Simulation Studies 

The simulation was used to investigate the closed-loop aircraft-autopilot problems 
and , subsequently , to define autopilot improvement modifications . The effects of 
static pressure angle-of-attack sensitivity and outer-loop gains on autopilot per- 
formance were the primary areas of investigation . 

For the study of the modified altitude hold autopilot , the SAS loop was not altered 
and autopilot hardware changes were kept to a minimum . In addition , because the 
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800 
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Figure 10. Original altitude hold operation with 
angle- 0 f-attack-insensitive static pressure, p^ . 


attitude hold mode of the autopilot performed well, the primary work on the autopilot 
was concentrated on the altitude control loops. 

The altitude hold investigation consisted of two goals: to improve the damping 
of the low frequency mode and to eliminate the high frequency mode . The first 
series of parametric studies concentrated on improving the characteristics of the low 
frequency mode with static pressure ports having no sensitivity to angle of attack. 
Once acceptable low frequency mode characteristics were obtained , the studies 
were devoted to determining a way to eliminate the high frequency mode. 
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Simulation excitation .— A number of methods were used to excite the aircraft- 
autopilot simulation . The primary methods were a step eleven input for a defined 
length of time; autopilot engagement from an initial rate of climb, which is typical 
of flight experience; and a 0.9° A0 step function input to the attitude loop of the 
autopilot , which is equivalent to a pilot pitch wheel input . The A0 step function 
was incorporated in the final hardware to provide an error signal of constant value 
to test autopilot response characteristics; therefore, the A0 step input was also 
used in the simulation for comparison purposes . 

In addition to the above disturbances, the autopilot's sensitivity to atmospheric 
pressure and temperature variations was investigated . 

Original altitude hold performance .— Because the lag and attenuation charac- 
teristics of the pitot- static system are not well defined, especially at the high 

frequencies of interest, the effective p (a) sensed by the ADC could vary from the 

s 

values shown in figure 5. Therefore, the simulation of p_(a) was varied by a factor 
until the data best matched the flight data, which occurred at a p„(a) of approxi- 

o 

mately one-half of the nominal setting . A time history of the simulation response of 
the original altitude hold autopilot to a 5-second, 2° elevon step with the ADC 
connected to p is presented in figure 11. 

Figures 12 and 13 are simulation responses of the original altitude hold mode 
with the ADC connected to p (no angle-of-attack sensitivity) . In figure 12, the 

disturbance (a 5-second, 2° elevon step) was small and the aircraft reached a 
neutral limit cycle with altitude variations of ±18 meters (±60 feet) . The limit cycle 
is due to the ADC threshold . 

In the response shown in figure 13 , the disturbance was an initial rate of 
climb of approximately 21 meters per second (70 feet per second) . The pitch SAS 
was engaged during the climb and , as the aircraft passed through an altitude of 
23,600 meters (77,500 feet) , the attitude and altitude hold autopilot modes were 
engaged . In addition , in the simulation runs with initial rates of climb , a constant 
2° excess power lever angle (PLA) from trimmed, Mach 3.0, 23,600-meter 
(77,500-foot) altitude flight conditions was used to partially minimize the Mach 
number loss . In this example , the disturbance was large and the aircraft entered 
a low frequency, saturated autopilot limit cycle. The 6^ simulation result was 

similar to the latter portion of the data in figure 10, which was obtained with the 
same autopilot-ADC configuration. 

Parametric autopilot gain variations .—The parametric studies presented in 
this section are perturbations about the selected gains with p (a) equal to zero. 

o 

All the gains are referred to in terms of a fraction of the nominal, or original, 
altitude hold autopilot gains . The disturbance was an initial rate of climb of 
approximately 21 meters per second (70 feet per second) with a 2° excess PLA and 
the autopilot engaged at an indicated altitude of 23,600 meters (77 ,500 feet) . 
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Figure 11. Simulation response of original altitude hold autopilot 
to a 5-second, 2° elevon step; p.(a) = 1/2 nominal. 

o 


The simulation responses were summarized in terms of first overshoot peak; 
period of the first oscillation; a damping index parameter, DI, defined in appendix B; 
and the time required to reach and maintain altitude within ±9 . 1 meters (±30 feet) 
of steady state, which is denoted as t^^. 

The gains investigated were , K„ , K„, and K. (fig. 6) . As each gain 

was varied , the remaining gains were held constant at the values selected for incor- 
poration in the modified autopilot . The final modified gain values in terms of a 
fraction of nominal were as follows: 
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Figure 12. Simulation response of original altitude hold mode 

to 5-second, 2° elevon step; p„(a) = 0. 
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K = 1/2 


K = 1 


^0 = 1 




The Kq parametric variation results are presented in figure 14. The first 

®/h 


overshoot is a constant 122 meters (400 feet) for all values of . The period 
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Autopilot 
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Time, min 

Figure 13. Simulation response of original altitude hold mode to 
initial rate of climb of approximately 21 m/sec (70 ft/sec); P^(a) = 0. 


is a constant 50 seconds for greater than one-fourth nominal; however, for low 

^fh 

values of the period becomes difficult to define because of the extremely long 

^fh 

time required to return to steady state . When K_ is zero , the altitude does not 

Vh 

return to the engage altitude, but converges to a constant bias. Both DI and 
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Figure 14. Variation of modified altitude hold mode charac- 
teristics with . In fractions of nominal settings, = 1, 

^fh 

Kq = 1/2, Kq = 1, and = 1. 
h 


reach their minimum values at the final modified gain; however, DI and degrade 
rapidly as K„ is increased or decreased . 

Vh 


The parametric variation results are presented in figure 15 . The first 

overshoot is a constant 122 meters (400 feet) for all values of K„ ; however, the 

period decreases from 90 seconds to 25 seconds as increases from one-fourth 

h 

nominal to nominal . Both DI and t^^ reach their minimum values at the final modified 

gain. Again, DI and degrade significantly as is either increased or 

60 


decreased . 
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Figure 15. Variation of modified altitude hold mode characteristics 
with K„ . In fractions of nominal settings, K„ = 1, - 1/4, 

Kq = 1, and Ky^ = 1. 


The Kq ^ parametric variation results are presented in figure 16. The first 
overshoot is 122 meters (400 feet) for all values of K„ ; however, the period 
increases from 40 seconds to 115 seconds as increases from one-half nominal 

to twice nominal . DI and t^^ are at or near their minimum values at the nominal 

gains . However , the degradation in characteristics is not as significant for small 
increases in as for the gains previously discussed . 
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Figure 16. Variation of modified altitude hold mode characteristics 
with . In fractions of nominal settings, Kr, = 1/2, = 1/4, 

^h ^fh 

Kq = 1, and = 1. 


The Kq gain variation in figure 17 illustrates the effect of the attitude loop on 

altitude hold . The parametric variation of Kq illustrates that the nominal attitude 

hold mode is a good inner loop for the modified altitude hold autopilot . The major 
effect of reducing Kq is to seriously degrade DI , whereas an increase in Kq 

significantly increases . 

The gain variation shown in figure 18 illustrates the effect of varying the 

modified , Ka > s^^d gains simultaneously . A reduction of K, increases 

Vh 

the first overshoot, the period, DI, and t^^; an increase in increases the period, 
DI, and t^^. 

Simulation response time histories of the modified autopilot configuration with 

p (a) equal to zero are presented in figures 19 and 20. Figure 19 is the response 
s 
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Figure 17. Variation of modified altitude hold mode characteristics 
with K„. In fractions of nominal settings, = 1, = 1/2, 
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Figure 18. Variation of modified altitude hold 
with In fractions of nominal settings, 


mode characteristics 
= = 1 , 


K. = 1/2, and - 1/4. 
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Figure 19. Simulation response of modified altitude hold to 5-second, 

2° elevon step; p (a) - 0. 

s 


to a 5-second, 2° elevon step, and figure 20 is an engagement of altitude hold from 
an initial rate of climb of 21 meters per second (70 feet per second) . In both 
cases , the low frequency mode is well damped , as was indicated in the parametric 
studies. (The flattening in the Ah trace in figure 20 and in the Ah and AM traces in 
other time histories represents the limit of the recording device . ) 

Effe ct of sensitivity of static pressure to angl e of attack.— Simulation time 

histories illustrating the effect of various factors of nominal p (a) on the original 

s 

autopilot are presented in figure 21. The autopilot was engaged from the start 

and the disturbance was a 5-second, 2° elevon step. As p (a) increased from zero 

s 

to nominal, the short-period mode converged to a large-amplitude limit cycle; 

as p„(a) decreased from zero to negative nominal, the long-period mode converged 
s 

to a large-amplitude limit cycle . 
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(a) p^(a) = 1. 



Time, min Time, min 


(b) p^(a) = 1/2. 



Time, min Time, min 


(c) Pgfa; = 0. 



Time, min Time, min 

Ps(a)=-l/2. 

Figure 21. Simulation responses of original altitude hold for various 

fractions of nominal p (a) . Input is 5-second, 2° eleven step. 
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Figure 21. Concluded. 


Simulation responses of the modified autopilot with p (a) values of one-half 

s 

nominal and nominal are presented in figures 22 and 23. The disturbance in both 

examples was a 5-second, 2° elevon step. The p (a) value of one-half nominal 

s 

(fig. 22) is representative of the majority of flight data; however, larger values 
of angle of attack (such as in a turn) would tend to drive p (a) toward the nominal 

o 

value (fig. 23) . A divergent-convergent low-amplitude oscillation exists for a 

p (a) of one-half nominal and, as p (a) approaches the nominal value, the oscil- 
s s 

lation gradually approaches a steady-state, large-amplitude limit cycle. 

Compensation for static pressure angle -of -attack sensitivity .—To compensate 
for the angle-of-attack sensitivity of static pressure , the high frequency oscilla- 
tions were damped by feeding back high-passed pitch rate . This method is rela- 
tively insensitive to nose boom configuration . 

The loop used to damp the high frequency mode consisted primarily of high- 
passed (t = 2) pitch rate summed with the attitude hold loop of the autopilot 
(fig. 24) . In addition, the high-passed pitch rate loop included a notch filter 
for aircraft bending compensation and a cutoff filter (i = 0.02) . The nominal 

was 1.07 deg 0/ (deg/sec 0) . Figure 25 is a simulation response to a 5-second, 

2° elevon step with the modified altitude hold gains and a nominal p (a) . As 

o 

shown, the high frequency mode of figure 23 is completely suppressed. The gain, 
Kx , of the high-passed pitch rate loop is one-half nominal in this example . Analyt- 

V7 

ically, larger values are desirable; however, larger values could feed significant 
system noise to the elevon . 

Impact pressure feedback .—At high altitude cruise conditions, the YF-12 air- 
plane has a phugoid mode with a damping range from slightly positive to negative . 
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Figure 22. Simulation response of modified altitude hold to 5-second, 

2° elevon step; p (a) - 1/2 nominal, 
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Figure 24. High-passed pitch rate loop used to compensate for 
static pressure angle-of-attack sensitivity (see fig. 6). 
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Figure 25. Simulation response of modified altitude hold to 5-second, 
2® elevon step. P^(<^) = nominal and = 1/2 nominal. 
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Because the original altitude hold was generally unsatisfactory , the pilot flew the 

constant altitude runs manually (SAS on but autopilot off) . One control concept to 

improve SAS -on handling qualities was to significantly increase phugoid damping 

by feeding back high-passed (t = 5) impact pressure, q . A filter (t = 1) was 

o 

2 

included in the feedback loop and the high-pass gain was -9.6 deg bj (N/m q^) 

(-0.2 deg 6 /(Ib/ft^ q )) . Figure 26 presents a comparison of the basic aircraft 

phugoid mode and the phugoid mode augmented by q feedback . The damping is 

o 

improved from a nearly neutral to a deadbeat condition . 


% feedback phugoid 



2000 


H 0 


Ah, 

ft 


-2000 



Figure 26. Simulation comparison of basic aircraft phugoid and q^- 
augmented phugoid modes for 10-second , 2^ elevon step; SAS on. 


The effect of adding the q loop to the modified altitude hold autopilot was 

c 

negligible, as shown by a comparison of figure 27 (q^ on) and figure 20 (q^ off) . 
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Figure 27. Simulation response of modified altitude hold to initial rate of 
climb of approximately 21 m/sec (70 ft/sec) with Pgfct) = 0 and feedback. 


Pitch Flight Control System Hardware Changes 

The pitch autopilot changes implemented based on the previous studies are 
presented in figure 28. The gyro used in the high-passed pitch rate loop was the 
backup for pitch SAS . An additional ADC was modified to provide . which 

is linear with q over the range of interest . In addition to the changes indicated 

on figure 28, the altitude error and integral altitude gains, K» and K_ , were 

reduced to one-half and one-fourth of their nominal values . The altitude rate 
gain , Kn , and Ka were left at their nominal values . The nominal value of K 

h c 
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From Mach and 
altitude loops 


Figure 28. Pitch autopilot modifications (see fig. 6). 



was -330 deg 6 /log-„ q . System flexibility was obtained by providing a gain 
%i X u o 

variation of at least ±6 decibels from the modified values for Kx , K„ , , Kq , 

® 

and Kq 
^j-h 


The cockpit flight test panel contained the on-off switch , the modified- 

original altitude hold select switch , and the five variable potentiometers . In the 
modified altitude hold position , the gains were switched from their original values 
to the modified values and the high-passed pitch rate loop was included . The A0 
disturbance function switch was also located on the flight test panel and could be 
engaged in any autopilot mode . 


In-Flight Performance of Modified Autopilot 

The first engagement of the modified altitude hold mode in flight introduced a 
hardware system problem that had not been anticipated . When the modified altitude 
hold mode was engaged , the system went into a limit cycle that produced an eleven 
output of ±1® at a frequency of 1 hertz (fig. 29) . The resultant normal accelerations 
of ±0 . 03g at the cockpit were evident to the crew . (Note that this is the only figure 
with cockpit acceleration; all other accelerations were measured at the center of 
gravity . ) 

Although the altitude was maintained accurately , the system performance was 
clearly undesirable . After approximately 2 minutes of the run , the pilot reduced 
the high-passed pitch rate gain, K^, from 1.07 to 0.65. An elevon limit cycle of 

±0.3° still persisted, but the crew could no longer distinguish the accelerations 
from the nominal level of aircraft vibration. was subsequently reduced to 0.50 

on this flight . Further attempts to reduce produced the following gain schedule 

as a function of Mach number: 




M < 1.4 . 

0.18 

1.4<M<2.75 . 

0.35 

M > 2.75 

0.50 


When the gain was reduced below the above values , the effects of the angle-of-attack 
sensitivity of static pressure became apparent . 

Appendix C analyzes the problem caused by the negative deficiency character- 
istic (ref. 10) of the backup pitch rate gyro. Briefly, the nonlinear gyro output 
around zero caused the rate signal to bounce from positive to negative in this 
region and effectively produced large gains at low input levels. Subsequent 
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Figure 29. Effect of pitch gyro nonlinearity (negative 
deficiency) in the high-passed pitch rate feedback loop. 

research showed that the gyro was noted for reliability rather than precision . It 
was decided that by using the gain schedule given previously , the crew would 

not be able to discern the limit cycle . Because of this and because the gyro did 
not affect overall autopilot performance , no attempt was made to improve or change 
the gyro. 

The following sections discuss the in-flight performance of the modified altitude 
hold autopilot , including the results of some limited changes in the modified control 
system . 

Modified altitude hold.— Typical flight test data obtained with the modified 
altitude hold autopilot are shown in figure 30. After engagement, the autopilot 
maintained altitude within ±8 meters (±25 feet) for the 4-minute duration of the run. 
The long-period (approximately 35-second) , low-amplitude oscillation is due to 
the threshold of the air data computer. After 2 1/2 minutes of the run, the pilot 


38 



II 


Autopilot 

Autopilot disenaaoed 



0 12 3 4 

Time, min 


Figure 30. Performance of modified autopilot. 

began a deceleration. As illustrated, there was no degradation of altitude hold, 
even though Mach number was decreased by 0.4. The high frequency (1-hertz) , 
low-amplitude (±0.02®) oscillation of the eleven is the short-period limit cycle 
described previously (fig. 29). 

A typical YF-12 experimental flight for purposes other than control system 
research consists of 1-minute periods of flight at stabilized Mach number and 
altitude conditions for as many conditions as can be scheduled in the flight plan . 
For these flights , pilots have found the modified autopilot valuable for rapidly 
attaining and maintaining altitude. Figure 31 illustrates a typical engagement 
of the modified autopilot with an initial rate of descent of approximately 400 meters 
per minute (1300 feet per minute) . 
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Figure 31. Performance of modified autopilot with 
initial rate of descent. 


The pilots described the atmosphere as stable for both modified autopilot runs; 
however , any latent aircraft- control system instabilities would probably have been 
excited by the deceleration in figure 30 or the initial rate of descent condition in 
figure 31 . 

A more exacting task for the modified autopilot is operation in turning flight , 
particularly the transition between wings level and steady-state turning flight. 

An example of a high speed turn with modified altitude hold is presented in fig- 
ure 32. The autopilot was engaged with wings level and, after 1 minute, a 
35° bank turn was entered. Overall, altitude control was within a 30-meter 
(100-foot) band. 
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Autopilot 



Time, min 

Figure 32. Modified altitude hold performance with turn 
entry. 35° bank angle. 

Figure 33 is a time history of wings level, turn entry, and steady-state turning 
flight . This figure illustrates a number of problem areas that are peculiar to the 
present YF-12 control system but must be considered in the design of any altitude 
hold autopilot. 

The altitude hold during the wings level portion was fairly good; however, the 
rapid turn entry excited a high frequency oscillation , which is caused by the static 
pressure angle-of-attack sensitivity. In this instance, the of 0.35 was lower 

than the Kx of 0.50 that was recommended for use at this Mach number . The auto- 

v7 

trim also reached its limit during this time period , thus directing the entire auto- 
pilot command to the pitch series servos through the ±2.3® limiter (fig. 6) . Because 
the low frequency autopilot command was not removed by the autotrim, the 2.3® auto- 
pilot authority limit was reached (after 6 1/2 minutes of the time history) and, at 
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Figure 33. Modified altitude hold performance with turn entry. Autopilot 
and autotrim saturation; 38° bank angle; - 0. 35. 


that point , altitude could no longer be maintained . An increase in autotrim 
authority would eliminate this problem . 

The effectiveness of the high-passed pitch rate feedback loop in suppressing 
the angle-of-attack sensitivity effect is illustrated in figure 34. In this example, 
modified altitude hold was engaged , was set at the desired value , and a A0 

step disturbance was used to excite the aircraft-autopilot system. This test, which 
was conducted for values of 0.35, 0.18, and 0, clearly showed the increased 

effectiveness of the high-passed pitch rate feedback with increased values of Kx . 

0 


42 




0 2 4 6 8 

Time, min 


Figure 34. Modified altitude hold performance for various values 
of K^; response to A0 disturbance function. 


Effect of K 


Q on modified altitude hold. 


-In the first few flights of modified 


altitude hold operation it was observed that the integrator path of the control 
system had a significant deadband, which occasionally was offset from zero. The 
deadband was a result of the integrator , which is mechanical , and had a significant 
breakout voltage , which in turn was aggravated by the reduction being 

made upstream of the integrator. As a result, the original breakout voltage, which 
was equivalent to an altitude error of ±5 to ±8 meters (±15 to ±25 feet) , was changed 
to an equivalent altitude error of ±18 to ±30 meters (±60 to ±100 feet) for a K_ 


of one-fourth the nominal value . 
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The initial attempt to reduce the breakout level was to increase K 


e 


Xh 


Although 


an increase in also decreases the damping , it was desirable to observe the 

®/h 

tradeoff. Figure 35 presents the results of doubling the gain from the modified 
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Figure 35. 


Modified altitude hold; = 1/2 nominal. 

^fh 


value to a of one-half the nominal value , which also effectively halves the 

^fh 

integrator breakout level. The steady-state altitude hold is very good; however, 
the apparent low damping associated with the large overshoot after the first minimum 
is undesirable and objectionable to the pilots. Basically, the pilots feel that the 
sign of the altitude error should not reverse during the initial capture period . This 
configuration was clearly less satisfactory than the first autopilot modification; 
therefore , K„ was reset to one-fourth the nominal value . 
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On subsequent flights , the integrator operation became erratic due to the inte- 
gration of stray voltages , thus making altitude control impossible . To salvage 
altitude control capability quickly , the integrator was grounded , which is equiva- 
lent to setting to zero . Previous simulator studies indicated that the damping 

Vh 

characteristics would be good but that any initial altitude rate at the time of engage- 
ment would prevent the steady-state altitude from returning to the original engage 
altitude. A typical example of modified altitude hold with set to zero (fig. 36) 

jh 

shows that the initial rate of climb at the time of engagement resulted in a steady- 
state altitude error of 61 meters (200 feet) . 


Autopilot Autopilot 



Time, min 


Figure 36. Modified altitude hold; = 0. 

^fh 
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was 


Integrator loop improvement .— Vihen the integrator was repaired, 

repositioned downstream of the integrator . In this position , the value of Kq had 

jh 

no effect on the threshold of the integrator , which remained fixed between ±5 and 
±8 meters (±15 and ±25 feet) . Flight evaluation of the system with the improved 
integrator loop showed good performance. The time history shown in figure 37 was 
obtained with the improved integrator loop and a Kg of one-fourth nominal . The 

^J'h 

altitude was quickly returned to its engage condition without an overshoot. 


Autopilot 



Time, min 

Figure 37. Modified altitude hold with improved integrator 
resolution; K„ = 1/4 nominal. 
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Altitude hold performance with alternate s tatic pressure source The modified 
altitude hold autopilot was also evaluated with the air data computer connected to 
the angle-of-attack-insensitive static pressure source, p . With this arrangement, 

®3 


the high-passed pitch rate loop was not required; therefore, was set to zero. 

These tests were also conducted during the period when the integrator was inop- 
erable and , therefore , K„ was zero . 

®/h 


Figure 38 presents a typical example of altitude hold in this configuration and, 
as would be expected, the performance is similar to that shown in figure 36 where 
the air data computer is connected to p „ and is zero . 
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Figure 38. 


Modified altitude hold with K 


Vh 


set to zero and ADC 


connected to p„ ; Ka = 0. 
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Effects of pressure altitude fluctuations on altitude hold. — One disadvantage of 
an altitude hold autopilot that uses air data as its only source of altitude information 
is that the control system chases pressure changes. An example of altitude hold in 
a fluctuating atmosphere is shown in figure 39. The pressure altitude fluctuations 
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Figure 39. Modified altitude hold response to indicated pressure 
altitude fluctuations . = 0; ADC connected to p^ ; = 0. 


are equivalent to ±9 meters (±30 feet) with a period of 23 seconds. The ride 
produced under these circumstances was unacceptable (a^ = ±0.25g) and as a 

result the system was disengaged. The large acceleration excursions are due to 

» 

the large indicated altitude rate signal , In p , provided to the autopilot by the 

o 

ADC in conjunction with the relatively large K_, on this signal . 
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Pressure fluctuations can create a severe problem when altitude hold control 
is used, as illustrated in figure 39. Although such occurrences are infrequent, 
the example illustrates problems that can be encountered and should be considered 
in the design stage . One way to minimize the problem would be to provide inertial 
information to the autopilot . A blended scheme could then use inertial information 
to control the high frequency component and air data information to control the low 
frequency component of autopilot control . 

Evaluation of impact pressure feedback . —Simulation studies showed that 

with SAS on, q feedback effectively damped the phugoid mode; however, with 
c 

altitude hold on , the effect of the q feedback was negligible . 

c 

The in-flight effectiveness of q feedback in damping the phugoid mode is 

c 

illustrated in figure 40. In this example, a 15-second thrust reduction was used 
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Figure 40. Effect of q^ feedback on phugoid mode. 

SAS on, autopilot off. 

49 



to excite the phugoid mode . The engagement of q was effective in damping the 

phugoid oscillation , which had a peak-to-peak altitude variation of nearly 
1220 meters (4000 feet) . 

Figures 41 and 42 present a comparison of altitude hold with and without the 
Qq feedback loop engaged . In both examples , the A6 disturbance function was put 

in and then taken out approximately 1 1/2 to 2 minutes later. Altitude hold is more 
sluggish with engaged and has a slight tendency to overshoot . A rougher ride 

is also produced by the q^ feedback , The rougher ride could be due to atmospheric 

temperature variations , which result in continual q variations . 
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Figure 41, Modified altitude hold with feedback engaged; 
response to A0 disturbance function. 
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Figure 42. Modified altitude hold without feedback; 
response to A0 disturbance function. 


Pilots' imp ressions of autopilot modifications .— The pilots' impressions of the 
modified altitude hold autopilot configurations ranged from very pleased to quite 
disappointed, depending on the control system configuration. The configuration 
most preferred by the pilots, based on all the autopilot experience, was the one 
that required the least amount of maneuvering from the engage point to steady 
state . If the altitude hold mode were engaged with the aircraft perfectly trimmed , 
nearly all the configurations would produce equivalent results . However , based 
on YF-12 flight experience, nearly all altitude hold engagements are from untrimmed 
conditions and have rates of climb . This initial rate of climb and the manner in 
which the control system reacts to it affect the pilot's opinion of the autopilot. 

The configuration most preferred by the pilots was that with the altitude integral 
gain set to zero (figs. 36 and 38) . An altitude error always exists when is 

®/h 
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zero and the system is engaged with a rate of climb . To the pilots , however , an 
altitude error is inconsequential as compared with extra maneuvering; when flying 
without the aid of the autopilot , the pilot accepts the altitude at which he initially 
stabilizes if it is within approximately 100 meters (300 feet) of the desired flight 
condition. Valuable time is wasted if a precise altitude is required. An altitude 
capture mode that would smoothly converge on a designated altitude would be 
desirable . 


The pilots' comments were favorable with the large integrator threshold, which 
resulted in well-damped control . The larger the threshold the better the pilots 
liked the system . 


Improving the integrator threshold (making it smaller) by moving the gain 
downstream (fig. 37) caused a degradation of pilot opinion and, in the case 

Vh 

where was increased , pilot comments were negative due to the low damping 

®/h 

characteristics (fig . 35) . 


Because the modified configuration with the improved integral loop threshold 
did not receive the most favorable pilot comments, was lowered to 0.15 of 

^fh 

the nominal value , at which point pilot comments indicated a definite improvement . 


The pilots were impressed with the ability of the feedback in damping the 
phugoid mode , although they detected no benefit in having feedback when altitude 
hold was engaged. 

Extensive handling qualities investigations were also made with and without 

q feedback for manual flightpath control (SAS on) and pitch autopilot control 
c 

(SAS on) . The effect of q^ feedback on handling qualities was negligible in all the 
configurations tested . 

With the altitude hold mode engaged and steady-state control established, the 
pilot's Mach hold task was easy as compared with simultaneous manual altitude and 
Mach control . 


Mach Hold Autopilot 

This section presents the results of the Mach hold autopilot investigation . The 
discussion covers the original Mach hold performance, the results of simulation 
studies with parametric variations , and the in-flight performance evaluation of the 
Mach hold mode with the ADC connected to the angle-of-attack-insensitive static 
pressure source. 
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Original Mach Hold Operation 


The Mach hold mode of the YF-12A autopilot was designed to operate over the 
entire Mach number range of the airplane; however, experience revealed long- 
period, large-amplitude divergent characteristics. Figure 43 presents a typical 
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Figure 43. Original Mach hold. 
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time history of a Mach hold engagement at approximately Mach 3 . 0 and an altitude 
of 23,600 meters (77,500 feet) with the air data computer connected to the angle- 
of-attack-sensitive static ports , p . Control of Mach number is good at first , 

with a period of 30 seconds and elevon inputs of ±1.5° commanding a normal acceler- 
ation of ±0 . 2g . Two minutes after engagement the oscillation has diverged suffi- 
ciently to saturate the autopilot output (A6^ = ±2.3°) , as evidenced by the ramp 

type of elevon control which is due to the autotrim function of the autopilot . The 
autotrim function begins to destabilize when the autopilot output saturates . In this 
configuration, a steady-state autopilot limit cycle develops from which the autopilot 
cannot recover with Mach hold engaged . The peak-to-peak altitude change was 
1370 meters (4500 feet) and the peak accelerations were greater than ±0.4g. There 
were no pilot inputs or obvious atmospheric disturbances to contribute to the 
problem . 


Simulation Studies 

The Mach hold autopilot problem was investigated briefly on the simulator . The 
areas of investigation were static pressure angle-of-attack sensitivity and outer- 
loop gain effectiveness . The primary means of disturbing the simulation was a drag 

2 2 

pulse, which consisted of a deceleration of 3.048 m/sec (10 ft/sec ) for a duration 
of 1 to 5 seconds and resulted in approximate Mach number losses of 0.01 to 0.04. 

A drag pulse is a controlled disturbance that can be duplicated on the aircraft by a 
thrust reduction and is also representative of a step temperature change . 

The Mach hold capability of the present control system is sensitive to variations 
of Pg(ct) > as shown in figure 44. For a 5-second drag pulse, which causes an initial 



(a) Pg(a) = 1. 

Figure 44. Simulation responses of original Mach hold for various 

2 

fractions of nominal p (a) , Input is a deceleration of 3. 048 m/sec 

s 

2 

(10 ft/sec ) for 5 seconds; and are nominal. 
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(b) Pg(a) = 3/4. 
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(c) p^(a) = 1/2. 
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(d) Pg(o.) = 1/4. 
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(e) Pg(a) = 0. 

Figure 44. Concluded. 


55 



Mach loss of 0.04, the system is unstable for a p (a) from three-fourths nominal to 

o 

nominal . As stated previously , the actual effective flight value of p (a) ranges 

from one-half nominal to nominal depending on angle of attack and attenuation at 
the higher frequencies. Elimination of p„(a) is obviously desirable. 

o 

In addition , the autopilot-aircraft nonlinearities make Mach hold stability 
sensitive to the disturbance magnitude. As shown in figure 45, the Mach hold 
control system with a nominal p (a) is stable for drag pulses of short duration . 

o 

However , for drag pulse durations of 4 seconds and longer , which result in Mach 
errors of 0.03 and greater, the system diverges. 
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(b) Drag duration = 2 seconds. 


Figure 45. Simulation responses of original Mach hold drag pulses 

2 

of various lengths . p (a) = nominal; deceleration = 3.048 m/sec 

s 

2 

(10 ft/sec ); Kq and are nominal. 
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(c) Drag duration = 3 seconds . 



(d) Drag duration = 4 seconds. 



(e) Drag duration = 5 seconds. 


Figure 45. Concluded. 


Simulation time histories for parametric variations of the outer-loop gains , 
Kq and , are presented in figure 46. For this study, p„(a) was set to 

one-half nominal to prevent the angle-of- attack sensitivity from obscuring the 
effects of the gain variations; the disturbance was a 5-second drag pulse. As 
shown , a reduction in K_ to one-half its nominal value minimizes autopilot 

^J'M 

saturation , thereby improving Mach control . As expected , an increase in 

fM 
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Figure 46. Simulation responses of original Mach hold for various 
outer-loop gain values, p (a) = 1/2 nominal; input is a deceleration 

o 

of 3. 048 m/sec^ (10 ft/sec^) for 5 seconds. 
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(d) K„ = 1 1/2 nominal, Kr, = nominal. 
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(e) Kr, = 1/2 nominal, Kr, = nominal. 

J'm 

Figure 46. Concluded. 

increases the divergence of the system. In addition, some increase in K^ appears 
to improve the stability of the system, whereas a decrease in K^ causes the system 
to diverge and saturate . 

A review of the simulation results emphasizes that at high speeds a conventional 
Mach hold autopilot with elevator control is sensitive to small Mach number changes , 
and that even normal atmospheric variations will cause a wavy ride, because Mach 
number is controlled through changes in altitude . 


Mach Hold Flight Results With ADC Connected to p 

As stated previously , the ADC was connected to the static pressure angle-of- 
attack-insensitive source , p , to eliminate the high frequency oscillations experi- 

3 

enced in the original altitude hold autopilot. Therefore, the influence of the p 

®3 

hookup on Mach hold was also investigated , although improvement in Mach hold was 
not a formal objective of the present study. 
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With the p hookup, the Mach hold mode was stable, as shown in figure 47, 

3 

and Mach number was held within ±0.02. However , to maintain Mach number , the 

Thrust 


Autopilot reduction 

engaged cycle 



0 12 3 4 

Time, min 


Figure 47. Mach hold operation with angle-of~attack~insensitive 
static pressure, p . 

aircraft was put into a climb (the power levers were not moved) because of the 
slight initial acceleration at the point of engagement . Approximately 2 minutes 
after Mach hold engagement, the pilot reduced thrust to produce a controlled 
deceleration disturbance function to demonstrate the Mach hold control system’s 
capabilities. As shown, the Mach number loss was minimal; however, a large 
altitude loss was necessary to compensate for the loss of thrust . In general , 
altitude changes required to control Mach number result in poor ride qualities . 

The autopilot signal was saturated for approximately 20 seconds during the 
recovery of Mach number , but the autopilot recovered and there was no evidence 
of long-period instabilities . 
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The gradual increase of altitude in the Mach hold example (fig. 47) emphasizes 
that a single controller does not permit simultaneous altitude and Mach control . 
Studies are underway for the implementation of an autothrottle Mach controller , 
which would permit the independent control of Mach number and altitude . 


CONCLUDING REMARKS 


An autopilot improvement program was undertaken to improve the poor altitude 
hold and Mach hold control system characteristics at high altitude, high speed 
flight conditions . 

Static pressure source sensitivity to angle of attack was found to have a signifi- 
cantly adverse effect on the altitude hold and Mach hold autopilot modes . 

Good altitude hold was demonstrated at high altitude using existing autopilot 
concepts and hardware with only minor modifications . These modifications consisted 
of adding a high-passed pitch rate feedback to compensate for the effects of angle- 
of-attack sensitivity on altitude measurements and reducing the altitude error gain 
and altitude integral gain. 

Reducing the integrator threshold and increasing the altitude integral gain 
improved the altitude tracking performance; however, these changes also caused 
unfavorable pilot comments due to an apparent sensed decrease in stability . 

Good Mach hold using elevator control was demonstrated at high Mach numbers 
with an angle-of-attack~insensitive static pressure hookup. However, large 
altitude changes were necessary to compensate for small Mach number changes , 
which resulted in poor ride qualities . 


Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif,, November 22, 1976 
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APPENDIX A. -DETAILED DESCRIPTION OF 


FLIGHT CONTROL SYSTEM ELEMENTS 


Nose Boom 

The air data information was obtained from a compensated nose boom as shown 
in figure 4 . The nose boom installation is described in detail in reference 7 . The 
p ports were connected to the air data computer and provided the primary altitude 

information to the pilot and the autopilot. The p ports provided pilot backup 

information . The p and p ports were located on the compensated portion of 

the nose boom to minimize the transonic static pressure position error correction. 

However , these compensated ports become sensitive to angle of attack as Mach 

number increases (ref. 11) . The Mach error, AM, can be converted to a static 

pressure error, Ap , by the expression 

s 

Ap = F(p /M)AM 
s s 

where F is a function of Mach number as illustrated in figure 48. The p ports 
were identical in location and configuration to the static ports on the standard NASA 



M 


Figure 48. Variation of Mach number. 
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APPENDIX A .-Continued 


pitot-static nose boom probe, which is described in detail in reference 12. Wind 
tunnel and flight data confirmed that the static pressure error due to angle of attack 
was negligible . A flight calibration of this probe position is presented in refer- 
ence 11 and, as would be expected, the transonic position correction for the p 

ports is much larger than that for the compensated ports, p and p 

®1 ®2 


The air data computer is approximately 6 meters (20 feet) aft of the total and 
static pressure ports and, therefore, a significant amount of lag is present in the 
indicated p. and p measurements. Ground calibrations of the p. system gave 

h ® ^2 


a sea level lag constant, , of 0.27 second, which can be converted to a 



time constant at a given Mach number and altitude by the expression 


P, "si 'Psl'Pt 

'2 Pf 


An in-flight determination of the p system lag constant was performed at 

®2 

the aircraft short-period frequency. The sea level lag constant, , was 



0.056 second, which converts to a time constant at the altitude of interest by the 
expression 


= T, 


(P07/P0 ) 


Air Da ta Computer 

The air data computer receives total and static pressure data from the compen- 
sated nose boom installation and computes Mach number and altitude information for 
pilot display and autopilot information . The ADC is an electromechanical device 
that uses two force-rebalancing pressure sensor servo loops . Figure 49 illustrates 

the operation of the static pressure, p , and impact pressure, q , servo loops, 

s c 

which are similar . A change of pressure within the transducer causes the evacuated 
bellows to exert a torque about a torsion bar, which, in turn, causes a servo loop 
to exert a rebalancing torque about the opposite end of the torsion bar . The twisting 
of the torsion bar is proportional to the pressure . The mechanism measures pres- 
sures with only minute bellows displacements, minimizing the adverse effects of 
nonlinear bellows spring rates and bellows hysteresis. In both loops, the pres- 
sure rates, p^ and q^, are measured by the velocity generator, which measures the 

rate of the motor displacement . Mach number , altitude , and associated rate infor- 
mation are determined by the use of cams . 
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APPENDIX A. -Continued 



Figure 49. Diagram of pressure computing servo loop. 


Flight data input and output characteristics of the ADC static pressure servo 
loop at Mach 3.0 and an altitude of 23,600 meters (77,500 feet) are presented in 
figure 50. These data are the result of a control system-induced limit cycle. The 
actual pressure variation input to the ADC is presented in figure 50(a) . The ADC 
output signals are voltages proportional to the altitude change from the engage 
point and to the altitude rate (fig. 50(b)) . The threshold is then determined by 
comparing the constant portions of either the altitude or altitude rate trace 
(fig. 50(b)) with the actual pressure altitude variation (fig. 50(a)) . At an 
altitude of 23,600 meters (77,500 feet) , the threshold is approximately 3.7 meters 
(12 feet) . 



Time, min 


(a) Air data computer input. 


Figure 50. Air data computer input and output characteristics . 
M«3.0, 23, 600 meters (77, 500 feet) . 
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(b) Air data computer outputs . 


Figure 50. Concluded. 


The frequency response characteristics of the ADC were determined in the lab 
and are presented in figure 51. The frequency response tests were perturbations 



u, rad/sec 


APc 


input 


N/m^ (ib/ft^) 
± 3.59 (± 0 . 075 ) 
± 7.18 (± 0 . 150 ) 
± 14.36 (± 0 . 300 ) 
± 23.94 (± 0 . 500 ) 
± 47.88 (± 1 . 000 ) 


Figure 51 . Frequency response of air data computer 
for various input levels. 
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APPENDIX A .-Continued 


about the nominal flight condition of Mach 3.0 and an altitude of 23,600 meters 
(77,500 feet) . A minimum amount of tubing was used to connect the pressure 
oscillator to the ADC . The frequency response is constant to approximately 
0.5 radian per second; however, the lower amplitude pressure variations are 
significantly attenuated by the threshold . 


Original Altitude and Mach Hold Control Syste ms 


The YF-12 pitch-axis vehicle control system is shown in figure 6. The time 
constant and gains , which are scheduled as functions of flight condition , are 
presented in figure 52 . 

0 200 400 600 800 1000 1200 1400 

1 ^ ^ ^ ^ ^ ^ I 

L’ 


degfig 
deg/sec e 

0 10 20 30 40 50 60 70 

kN/m^ 

C’ 

(a) Schedule 1: pitch rate -to -eleven gain schedule. 



logio Pj (Pj in in. Hg) 

1.6 1.2 .8 .4 0 -.4 -.8 



iogioPj <Pj in N/m^) 


(b) Schedule 2: pitch attitude low-pass time constant. 


Figure 52. Autopilot schedules. 
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logic ''c 


(c) Schedule 3: pitch attitude gain schedule. 


logjoPs (Pj in in. Hg) 

.8 .6 .4 .2 0 -.2 -.4 -.6 



logic '’s 

(d) Schedule 4: altitude displacement gain schedule. 


Figure 52. Continued. 



APPENDIX A. -Continued 

logjoPj <P; in in. Hg) 

■ 8 .6 .4 .2 0 -.2 -.4 -.6 



logjO Pj (Ps in N/m^) 


fej Schedule 5: altitude integral gain schedule. 


30 


deg/sec 0 
Aln(R-l) 

0 

5.4 5.0 4.6 4.2 3.8 3.4 3.0 2.6 

logio Pj (P5 in N/m^) 

(f) Schedule 6: Mach integral gain schedule. 


logioPj (Pj in in. Hg) 
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I ^ ^ ^ I ^ I 



logioPj (Pj in in. Hg) 

1.6 1.2 .8 .4 0 -.4 -.8 
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(g) Schedule 7: Mach displacement gain schedule. 


Figure 52. Continued . 
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(h) Schedule 8: Mach trim system compensation. 
Figure 52. Concluded. 


Pitch SAS is used full time on the aircraft and consists of proportional pitch 
rate feedback for short-period damping augmentation. At altitudes above 
15,200 meters (50,000 feet) , lagged pitch rate is added to increase the aircraft's 
effective static stability . In the SAS mode , speed stability augmentation is provided 
from Mach 0 . 2 to Mach 1.5. 

The primary autopilot mode is attitude hold; the altitude and Mach hold modes 
are outer loops of the attitude hold mode . Attitude hold consists of lagged pitch 
attitude and high-passed pitch attitude , which are summed . The signal is multiplied 
by a scheduled gain and then simultaneously sent to the pitch series servo and the 
autotrim and its compensator. There is an eleven limit of ±2.3° on the SAS servo 
input . The threshold of the autotrim is 0 . 1 degree per second . The pilot can 
command pitch attitude changes through a pitch attitude trim wheel . 

The altitude hold mode receives altitude and altitude rate information from the 
ADC in terms of Ain p and In p respectively. The quantity Ain p is a pertur- 

S S o 

bation around the flight condition that exists at the time of engagement . This 
perturbation is obtained by means of a clutch system within the ADC . The range of 

Ain p is equivalent to approximately ±155 meters (±507 feet) . At an altitude of 

s 

23,600 meters (77,500 feet) , the natural log quantities given above are related to 
altitude by the following expressions . In SI Units , 

Ah = (-6490)Aln p 

s 

It. = (-6490)ln p 

s 

In U . S . Customary Units , 

Ah = (-21,300)Aln p 

o 

fi = (-21,300)Tn^ 
s 
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APPENDIX A . -Concluded 


Altitude rate , altitude error , and integral altitude are then summed and added to 
the attitude loop . 

The Mach hold mode receives Mach information from the ADC in terms of 
Ain (R - 1) , which is a perturbation around the flight condition at the time of 
engagement and a clutched function of the ADC . The range of the AZn (R - 1) is 
approximately equivalent to ±0,25 Mach number at Mach 3.0. At Mach 3.0, 

Ain iR - 1) is related to AM by the expression 

AM = (1.43)AZn (i? - 1) 

Mach error and integral Mach error are then summed and added to the attitude 
loop , Mach rate feedback was originally designed into the system , but the signal 
was found to have excessive noise and did not contribute to Mach stability . 
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APPENDIX B.— DEFINITION OF DAMPING INDEX AND RELATIONSHIP 
WITH SECOND-ORDER SYSTEM SUBSIDENCE RATIO 


The subsidence ratio for a second-order system is defined as 



following sketch (ref. 13) . 


based on the 



The quantity X(t) can be expressed as 

= A 


X(t) 


cos 


(V) 


where A is an arbitrary amplitude. The time, t, at the peaks of the oscillation is 
expressed as 


1 A. ^ 

= |-;r + I 

03 A2 / 


Evaluating X(t) at peak values of the oscillation yields the following subsidence 
ratio definition: 



&%m 

e 


where 8 = ^ • 

The subsidence ratio definition assumes that the parameter measured is linear 
and oscillating about zero or some known bias , which must be removed prior to 
evaluation. The simulation used in the autopilot studies was nonlinear and greater 
than second order; therefore , a gpven parameter did not oscillate about a known 
bias . However , much of the basic response information approximated second-order 
data superimposed on an exponential curve . 
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APPENDIX B.- Concluded 


Based on the following sketch, a damping parameter similar to the subsidence 
ratio was defined . This damping parameter , referred to herein as the damping 
index, DI, is given by the following expression: 


DI = 




When this formula is applied to a second-order system , 


which is identical to the subsidence ratio definition when m equals 1 . 
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APPENDIX C.— ANALYSIS OF MODIFIED AUTOPILOT 


HIGH FREQUENCY OSCILLATION 


Experimental studies showed that measured static pressure was sensitive to 
angle of attack and that this sensitivity tended to reduce the short-period stability 
with the outer loops engaged . To partially compensate for this destabilizing effect , 
high-passed pitch rate was implemented in the modified altitude hold autopilot as 
an additional feedback loop . 

Figure 53 is a block diagram of the attitude portion of the modified system . Pitch 
rate is sensed by a backup SAS gyro , which is forward of the center of gravity; 

50s 

high-passed through the filter (g + q 5 ) (gi : 50 ) ’ then summed with the autopilot 
outer-loop signals at the gain schedule 3 input summer (fig. 28). 


High-pass loop 



Figure 53. Block diagram of modified autopilot attitude hold mode. 
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APPENDIX C. -Continued 


On the first flight after implementation of the modified autopilot , with the 
high-passed pitch rate feedback loop engaged and set at a gain of 1.07, the 
aircraft developed an undesirable residual oscillation . The oscillation , which 
had a frequency of approximately 1 hertz , was bounded at a peak-to-peak pitch 
rate of less than 0 . 6 degree per second . 

The transfer functions that describe the pertinent elements of the system modeled 
(fig. 53) are as follows: 


B , . _ 

o o O 

« 4- cOg 


where 15.7 < < 20.6, = 0.016, and = 0.05; and 


0 ^ ^ 
^(s) = 


0.575s(s + 0.008) (s + 0.1646) (s - 0.0 049) 

(s - 0.00278) (s + 0.0000356 ± 0.055j)^(s + 0.133 ± 1.20/)^ 


1470 


" 72 


AAs) = 

s + 46s + 1470 
A(N) = ±0.15° hysteresis 
2210 


(s + 75.4s + 2210) (1 + 0.05s) 

F (s) = 0175 

(s + 0.175) 


F (s) = — O; ^0^ 

^21^1 (s + 0.4) 

F (s) = — — — — ^0^ 

HP'' 1 (s + 0.5)(s + 50) 


, 0- ^ 4S . 400) 

(s + 40s + 400) 

GyiN) = ±0.1 deg/sec negative deficiency 


- 0-15(s + 8) , 4(s + l) 

^ (s + 4) (s + 10.4) (s + 0.077) 


8„ =3.6 


74 



APPENDIX C .-Continued 


Omitting the transmission linkage nonlinearity , A (N) , the system was manipu- 
lated to an equivalent form, which is presented in figure 54. In the equivalent 
system, G' is) combines the SAS and attitude closed-loop transfer functions. To 
isolate gyro nonlinearity , high-passed pitch rate is depicted as the only feedback 
loop. In addition, the filter zeros of F^^(s) are assumed to cancel the bending mode 

characteristics . 


High-passed pitch rate loop 



®c G'(s) 

e " 1 + K.G'(s)H'(s)G^(N) 

10 Y 


Figure 54. Equivalent modified attitude hold system. 


A linear analysis of the equivalent system (fig. 54) was performed with Gy(N) 
assumed to be 1. The open-loop frequency response, K^G' is)H' (s) , is presented 
in figure 55 . The system shows an increase in 9/8^ of 40 decibels per decade at 



Figure 55. Open-loop frequency response of linear modified attitude 
hold mode including high-passed pitch rate feedback loop. 
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APPENDIX C. -Continued 


the lower frequencies and attenuation of approximately 160 decibels per decade at 
the higher frequencies . In addition , the expected phase shift accompanies the 
amplitude variation. For the approximate design gain setting, = 1.0, in the 

high-pass loop, the system has a phase margin of 45° and a gain margin of approx- 
imately 10 decibels. In general, the predicted open-loop margins are reasonable 
for this type of system. Additional gain or phase lag, or both, is required to 
sustain limit cycles . After considering possible sources , such as transmission 
linkage hysteresis , gyro location , and gyro nonlinearity , the gyro was believed to 
be the element most likely to provide the necessary gain and phase lag . 

The gyro was removed from the aircraft and mounted on a rate table , and normal 
oscillations and rate tests were conducted. The input to the rate table was at a 
frequency of 1 hertz and amplitudes slightly above those normally experienced in 
flight. As shown in figure 56, when the rate is approximately zero, the output is 



Figure 56. Negative deficiency of gyro used in high- 
passed pitch rate feedback loop, f = 1 hertz. 


particularly different and noticeable as compared to the rest of the cycle . This type 
of nonlinearity , referred to as negative deficiency , is common to rate sensors , 
occurring normally with increasing age. Presented in figure 57 are the results of 



o Increasing 
□ Decreasing 


Gyro input, deg/sec 


Figure 57. Negative deficiency characteristic of gyro. 
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APPENDIX C , —Continued 


constant-rate tests . The maximum input rates are typical of those measured in 
flight . The input to the rate table is a complete continuous cycle , first increasing 
in the clockwise direction , then decreasing and increasing through the counterclock- 
wise direction , and finally returning to zero . From the constant-rate tests it is 
evident that conservative values for the amount of deadband would be within 
±0 . 1 degree per second at low rates . At the higher rates , above 1 degree per 
second, the sensed rate of 1.7 mV/deg/ sec shows good agreement with the design 
scale factor. The increase in gain and phase lag descriptive of this type of non- 
linearity is presented as a describing function in figure 58 , which is adapted from 



0 |- 


-10 


-20 - 
<P. 
deg 

-30 - 


-40 


-50 >- 



Figure 58. Negative deficiency describing function (adapted from ref. 10) . 


reference 10. Classical system stability analysis was used with the describing 
function to predict the limit cycle frequency and amplitude . For example , with 
reference to figure 54 , the nonlinear gyro is manipulated to be cascaded with the 
other elements into a single loop . Direct application of the Nyquist stability criterion 
then provides insight regarding the cause and degree of this type of problem . 


The Nyquist criterion is essentially concerned with the condition when 
K^G' (s)H' (slGyCN) equals -1 or when the characteristic equation equals zero. 

Therefore, when this condition is just satisfied (the magnitude equals 1 and phase 
equals -180°), the system will be critically stable and able to sustain limit cycle 

oscillations . The two functions , K^G' (s)H' (s) and ^ , are separated for 
determining a solution. The open-loop transfer function, K^G' (s)H' (s) , is pre- 
sented in figure 55 as a function of frequency. The function ^ negative 

reciprocal of the describing function presented in figure 58 . In figure 59 , the two 
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u. 



(p, deg 


Figure 59. Phase versus amplitude of aircraft with modified 
altitude hold mode and nonlinearity (negative deficiency) 
in high-passed pitch rate feedback loop . 

functions are plotted as amplitude versus phase and are varied parametrically 
with frequency and amplitude , respectively . The intersection of the two curves , 
which occurs at a frequency of approximately 7 radians per second and a pitch rate 
of approximately 0.15 degree per second, predicts the frequency and amplitude 
of the limit cycle for a of 1 . 0 . It is evident that if the gain were reduced 

sufficiently , no intersection would exist; consequently , no residual oscillations 
would be predicted . 


Subsequent flight tests and simulator studies were conducted at the same flight 
test conditions. Figure 60 is a summary and comparison of data from the flight 



K^, deg/sec 

Figure 60. Modified autopilot attitude hold mode peak-to-peak 
oscillations and frequency versus high-passed pitch rate gain. 
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APPENDIX C.- Concluded 


tests , simulator studies , and limit cycle calculations . The calculations and simu- 
lator studies predicted lower peak-to-peak amplitudes than those experienced in 
the flight tests. These differences are believed to be caused by the ±0. 15° trans- 
mission linkage hysteresis between the inboard and outboard elevons . This non- 
linearity was not accounted for in the simulator studies or in the limit cycle 
calculations . 
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